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STABILITY OF STEADY STATES FOR A CLASS OF FLOCCULATION EQUATIONS 

WITH GROWTH AND REMOVAL 

INOM MIRZAEV* AND DAVID M. BORTZ*t 

Abstract. Flocculation is the process whereby particles (i.e., floes) in suspension reversibly combine and separate. The 
process is widespread in soft matter and aerosol physics as well as environmental science and engineering. We consider a general 
size-structured flocculation model, which describes the evolution of floes in an aqueous environment. Our work provides a 
unified treatment for many size-structured models in the environmental, industrial, medical, and marine engineering literature. 
In particular, our model accounts for basic biological phenomena in a population of microorganisms including growth, death, 
sedimentation, predation, renewal, fragmentation and aggregation. Our central goal in this paper is to rigorously investigate 
the long-term behavior of this generalized flocculation model. Using results from fixed point theory we derive conditions for the 
existence of continuous, non-trivial stationary solutions. We further apply the principle of linearized stability and semigroup 
compactness arguments to provide sufficient conditions for local exponential stability of stationary solutions as well as sufficient 
conditions for instability. 

The end results of this analytical development are relatively simple inequality-criteria which thus allows for the rapid 
evaluation of the existence and stability of a non-trivial stationary solution. To our knowledge, this work is the first to derive 
precise existence and stability criteria for such a generalized model. Lastly, we also provide an illustrating application of this 
criteria to several flocculation models. 
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populations dynamics, semigroup theory 

AMS subject classifications. 35Q02, 35P02, 45C02, 45G02, 45K02, 92B05 


1. Introduction. Flocculation is the process whereby particles (i.e., floes) in suspension reversibly 
combine and separate. The process is widespread in soft matter and aerosol physics as well as environmental 
science and engineering. A popular model for flocculation is a ID nonlinear partial integro-differential 
equation which describes the time-evolution of the particle size number density. In the engineering literature, 
this equation can be derived using a so-called population balance equation (PBE) framework and we direct 
the interested reader to the book by Ramkrishna [5DJ for more on this framework. 

Previous analytical work on these models focused on classes of flocculation equations that did not allow 
for the vital dynamics (i.e., birth and death) of individual particles. These phenomena are obviously critical 
features in the modeling of microbial flocculation. Accordingly, in this work we consider a general size- 
structured flocculation model which accounts for growth, aggregation, fragmentation, surface erosion and 
sedimentation. The variable p(t, x) denotes the number density of floes of size x at time t, and for a given 
interval U C (a.’o, aq), the function xu represents the characteristic function of the interval U. A floe is 
assumed to have a minimum and maximum size Xq and aq such that 0 < Xq < aq < oo.The equations for 
the microbial flocculation model can be written as 


where 


Q denotes growth 


d t p = F\p\ 

F\p\ := Q\p) + A[p\ + B\p], 

g\p] ■= - d x (gp ) - fl{x)p(t, x ), 


( 1 . 1 ) 


( 1 . 2 ) 


A denotes aggregation 


^ rx—xo 

A\p\ := -zX[ 2 x 0 ,x 1 ){x) / k a (x - y , y)p(t , x - y)p(t , y ) dy 

^ Jx o 

rxi — x 

- X[x 0 ,x 1 -x 0 ]P(t, x) / k a (x, y)p(t, y) dy 


(1.3) 
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and B denotes breakage 


r Xi l 

B\p] '■= X(xo,xi-xo]W J y)kf{y)p{t, y)dy - -X[ 2 x 0 ,x 1 )kf{x)b{t, x). (1.4) 

The boundary condition is traditionally defined at the smallest size xq and the intial condition is defined at 
t = 0 


g{x 0 )p(t, x 0 ) = / q(x)p(t, x)dx, p( 0 , x) = p 0 (x) e L 1 ^). 


Note that in vivo, there are no floes of size zero and the floes cannot grow indefinitely, so the only 
biologically realistic case is 0 < Xq < X\ < oo. However, when xq > 0, the characteristic functions appearing 


in equations (1.31 and (1.41 make our theoretical development rather cumbersome. Hence, for the sake of 


convenience, in this paper we consider the case xo = 0 < X\ < oo, and postpone the analysis of the case with 
0 < xq < X\ = oo for our future papers. Therefore, we will denote the closed interval by I := [0, X\\ and 
will make extensive use of this interval in our development. We carry out the analysis of this work (unless 
otherwise specified) on the space of absolutely integrable functions on /, denoted by T 1 (/)- We also note 
that well-posedness of the flocculation model on this space has been established by Banasiak and Lamb [S] . 

1.1. Background and model terms. In this section, we will provide a brief background and overiew 
of the individual terms in the general flocculation equation above. 


To begin, the Sinko-Streifer m terms in ( | 1 . 2 | ) correspond to the growth and removal of floes, respectively. 
The function g(x) represents the average growth rate of the floes of size x due to mitosis, and the coefficient 
g(x) represents a size-dependent removal rate due to gravitational sedimentation and cell death. Specifically, 
when an individual cell in the floe of size x divides into daughter cells, the new cells can remain with the floe, 
contributing in a increase in its total size. Conversely, a daughter cell can also leave the floe to form a new 
single-cell floe. This second case is modeled by McKendrick-von Foerster type renewal boundary conditions, 


g(x 0 )p(t, x 0 ) = / q(x)p{t, x)dx , 


(1.5) 


where the renewal rate q{x) represents the number of new cells that leave a floe of size x and enter the single 
cell population. We note that this boundary condition could also be used to model the surface erosion of floes, 
where single cells are eroded off the floe and enter single cell population. The well-posedness and stability 
of equilibrium solutions of the Sinko-Streifer equations has been extensively studied by many researchers 
using a wide variety of mathematical conditions [I55[ 134) HOI '40, 451 [15] . For numerical simulation of the 
model, a convergent numerical scheme has been proposed in ED3, and inverse problems for estimation of the 
parameters of the model have been discussed in mmm- 

The aggregation of floes into larger ones is modeled in by the Smoluchowski coagulation equation. 

The function k a {x, y) is the aggregation kernel, which describes the rate with which the floes of size x 
and y agglomerate to form a floe of size x + y. This equation has been widely used, e.g., to model the 
formation of clouds and smog in meteorology m, the kinetics of polymerization in biochemistry ESI, the 
clustering of planets, stars and galaxies in astrophysics m, and even schooling of fish in marine sciences j 4*5>j . 
The equation has also been the focus of considerable mathematical analysis. For the aggregation kernels 
satisfying the inequality k a {x, y) < 1+x+y, existence of mass conserving global in time solutions were proven 
mmm (f° r some suitable initial data), conversely, for aggregation kernels satisfying ( xy ) 7 / 2 < k a {x, y) 
with 1 < 7 < 2 , it has been shown that the total mass of the system blows up in a finite time (referred 
as a gelation time) [25] . For a review of further mathematical results, we refer readers to review articles 
by Aldous [5], Menon and Pego [43], and Wattis [541 and the book by Dubovskii [20]. Lastly, although the 
Smoluchowski equation has received substantial theoretical work, the derivation of analytical solutions for 
many realistic aggregation kernels has proven elusive. Towards this end, many discretization schemes for 
numerical simulations of the Smoluchowski equations have been proposed, and we refer interested readers to 
the review by Bortz m § 6 ]. 


The breakage of floes due to fragmentation is modeled by the terms in (1.41, where the fragmentation 
kernel kf{x) calculates the rate with which a floe of size x fragments. The integrable function T(x',y) 
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represents the post-fragmentation probability density of daughter floes for the fragmentation of the parent 
floes of size y. The post-fragmentation probability density function T is one of the least well-understood 
terms in the flocculation model. Many different forms are used in the literature, among which normal and 
log-normal densities are the most common |35j [53] . Recent modeling and computational work suggests that 
normal and log-normal forms for T are not correct and that a form closer to an arcsin(:r; y) density would 
be more accurate nsnni. However, in this work we do not restrict ourselves to any particular form of T, 
and instead simply assume that the function T satisfies the mass conservation requirement. In other words, 
all the fractions of daughter floes formed upon the fragmentation of a parent floe sum to unity, 


ry 

/ T(x; y) dx = 1 for all y £ (xq, x\\. 

J Xn 


1.2. Overview of model assumptions. The flocculation model, presented in (1.1 1 , is a generalization 


of many mathematical models appearing in the size-structured population modeling literature and has broad 
applications in environmental, industrial, medical, and marine sciences. For example, when the fragmentation 
kernel is omitted, kf = 0 , the flocculation model reduces to algal aggregation model used to describe evolution 
of phytoplankton community [3j. When the removal and renewal rates are set to zero, the flocculation model 
simplifies to a model used to describe the proliferation of Klebsiella pneumoniae in a bloodstream m- 
Furthermore, the flocculation model, with only growth and fragmentation terms, was used to investigate the 
elongation of prion polymers in infected cells ED US]. 

The flocculation model in this form (1.1 1 was first considered by Banasiak and Lamb in [2., where they 
employed the flocculation model to describe the dynamical behavior of phytoplankton cells. The authors 
showed that under some conditions the flocculation model is well-posed, i.e., there exist a unique, global in 
time, positive solution for every absolutely integrable initial distribution. For the case Xi = oo, Banasiak 
[§] establishes that for certain range of parameters, the solutions of the flocculation model blow up in finite 
time. Nevertheless, to the best of our knowledge, for the case x\ < oo the long-term behavior of this model 
has not been considered. This is mainly due to nonlinear nature of Smoluchowski coagulation equations 
used for modeling aggregation. Hence, our main goal in this paper is to study the long-term behavior of the 


broad class of flocculation models described in (1.1 1 . For the remainder of this work , we make the following 
assumptions 


(Al) 

9&C\l) 

9(x) 

> 0 for x £ I 

(A2) 

k a £ W ^“(J 

x /), 

k a {x, y) = k a (y, x) 


and k a (x, y) 

= 0 

if x + y > x\, 

(A3) 

9 € C(I) 

and /i 

> 0 a.e. on I , 

(A4) 

q £ l°°(i) 

and 

q > 0 a.e. on I , 

(A5) 

k f £ C{I) 

M o: 

1 = 0 and kf > 0 a.e. on I , 

(A 6 ) 

F(-, y) £ W 1 ’ 

00 oo, 

r(^; y) > 0 for x £ ( 0 , y] 


and r(x; y ) = 

= 0 for 

x £ (y, xi). 


Assumption (Al) states that the floe of any size has strictly positive growth rate. This in turn implies 
that floes can grow beyond the maximal size aq, i.e., the model ignores what happens beyond the maximal 
size aq (as many authors in the literature have done HU 02 [□ [26]). We also note that the Assumption (Al) 
generates biologically unrealistic condition g{ 0 ) > 0 , i.e., the floes of size zero also have positive growth rate. 
However, this assumption is crucial for our work, and thus we postpone the analysis of the case g(0) = 0 for 
our future papers. Assumption (A2) states that for the aggregates of size x and y the aggregation rate is 
zero if the combined size of the aggregates is larger than the maximal size. Lastly, Assumption (A3) on p(x) 
enforces continuous dependence of the removal on the size of a floe and ensures that every floe is removed 
with a non-negative rate. 

When the long-term behavior of biological populations is considered, many populations converge to 
a stable time-independent state. Thus, identifying conditions under which a population converges to a 
stationary state is one of the most important applications of mathematical population modeling. It is trivially 
true that a zero stationary solution exists, but we are also interested in non-trivial stationary solutions of 
the flocculation model. Hence, in Section [2] we first show that under some suitable conditions on the model 
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parameters the flocculation equation has at least one non-trivial (non-zero and non-negative) stationary 
solution. Once a stationary solution to a model is shown to exist, the next natural question is whether it is 
stable or unstable. When the associated evolution equation of a population model is linear, many of stability 
properties can be deduced from the spectral properties of this linear operator |18l [32] . However, almost no 
information about the operator can be deduced from the spectrum of a nonlinear operator |7j. Moreover, there 
is no general consensus among mathematicians on how to define spectrum of a nonlinear operator. Thus, 
our stability analysis in this work is based on the principle of linearized stability for nonlinear evolution 
equations [55] E3- Hence, in Section [3] we summarize the principle of linearized stability and linearize the 
flocculation model around its stationary solutions. In Section[3]we also derive conditions for the regularity of 
the linearized flocculation model. Next, in Sections [4] and [5] we derive sufficiency conditions for the linearized 
stability and instability of zero and non-zero stationary solutions. In Section [ 6 ] we illustrate our results with 
several examples. Finally, in Section [7] we summarize and discuss the conclusions of this work. 


2. Existence of a positive stationary solution. The flocculation model under our consideration 
( 1 . 1 ), accounts for physical mechanisms such as growth, removal, fragmentation, aggregation and renewal 
of microbial floes. Thus, under some conditions, which balance these mechanisms, one could reasonably 
expect that the model possesses a non-trivial stationary solution. Hence, our main goal in this section is to 
derive sufficient conditions for the model terms such that the equation ( 1 . 1 ) engenders a positive stationary 
solution. 

Recall that at a steady state we should have 


Pt = 0 = X\p\ ■ 


( 2 . 1 ) 


By Assumption (Al), we know that 1/g £ C(I ) and thus we can define p = f/g for some f £ C(I). The 
substitution of this / into ( 2 . 11 , integration between 0 and an arbitrary x, and rearrangment of the terms 
yields 


f(x) = [ 
Jo 




g(y) 


+ 


l 


r kg (z - y, y) 

I o g(z-y)g{y) 


f{z-y)f{y)dydz - f 
Jo 


JO Jz 

M 

g(z) 


1 r(z; y)kf(y) 

g{y) 

rXl k a (z, y) 

g(y) 


f(y) dy dz 


f(y) dydz. 


r 0 


We now define the operator <E> as 


$[/](x) := 


r " q df\miy-r k -MllpM mdy 

3 9\y) Jo 


T(z; y)k f (y) 


g{y) 


r kg{z - y, y) 
I o g{z-y)g(y ) 


f(z~y)f{y)dydz~ [ 
Jo 


Jo j ; 

m 

g{z) 


g(y) 

r xi k a (z , y) 

g(y) 


f(y) dydz 


f(y) dy dz. (2.2) 


'o 


and will use a fixed point theorem to prove the existence of a fixed point / of $ . This in turn will allow us 


to claim that equation ( 2.11 has at least one non-trivial positive solution. 


The use of fixed point theorems for showing existence of non-trivial stationary solutions is not new in size- 
structured population modeling. For example, fixed point theorems, based on Leray-Schauder degree theory, 
have been used to find stationary solutions of linear Sinko-Streifer type equations [75] HB1. Moreover, the 
Schauder fixed point theorem has been used to establish the existence of steady state solutions of nonlinear 
coagulation-fragmentation equations ]551 . For our purposes we will use the following fixed point theorem, 
and refer readers to [6] for the full discussion of the proofQ 

Theorem 2.1. Let X be a Banach space, K C X a closed convex cone, K r = K D B r (0), $ : K r —> K 
continuous such that $ ( K r ) is relatively compact. Assume that 

1. $[x] Ax for all ||x|| = r and A > 1. 

2. There exists a p £ { 0, r) and k £ AT\{0} such that 

x — $[x] ^ A k for all ||x|| = p and A > 0. 


1 Hereafter, we refer to the following theorem as “the fixed point theorem” 
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Then $ has at least one fixed point Xq £ K such that p < ||a;o|| < r. 


Next, we show that the operator $ defined in (2.2) satisfies the assumptions of the above theorem, which 


in turn implies existence of a positive stationary solution of the operator T. Since we have been working 
on the space of absolutely integrable functions on /, a natural candidate for the Banach space X would 
be However, to obtain sufficient regularity for stability analysis of a stationary solution, we choose 

X = C(I) with usual uniform norm ||-|| on I. We also denote the usual essential supremum of a function 
by ||■ Hoo- Since the positive cone in C(I), denoted by (C(I)) + , is closed and convex, we choose K to be 
(C(/)) + . Then K r = KnB r ( 0), where B r ( 0) C X is an open ball of radius r and centered at zero, and r has 
yet to be chosen. We are now in a position to state the main result of this section in the following theorem. 

Theorem 2.2. Assume that the conditions 


0 < q(x) + \k f (x) - p(x), 


(Cl) 


and 




T(*; y)dz + q{y) 



1 

9 (y) 


g(y) + 2 k f(y) ~ v(.y) 


dy < i, 


(C2) 


hold true for all x £ I. Then the operator <f> defined in (2.2 1 ) has at least one non-zero fixed point, /* £ K 
satisfying 


0<V< II/* IL < r 


(2.3) 


for some 77 , r > 0. Moreover, the non-zero and non-negative function 

P* = j£ <?(/) 


(2.4) 


is the stationary solution of the flocculation model defined in [1.1). 
Proof. For / £ K r we have 


$[p] > 


+ 


X1 qty± 

g(y) 

pX pX 

1 0 j z 

f(z) 


f(y) dy — / 

JO 

r(~; v)kf{y) 

g(y ) 


: kf{y)/2 + p(y) 

y(y) 

f(y) dy dz — [ 


f(y ) dy 

m r 


g(z) 




> 


g(z) 

m 

g{z) 


kg(z, y) 

g(y ) 


kq{z, y) 

g{y) 

f(y) dy 


f{y) dy dz 


dz + 


kf(y)f(y) 

g(y) 


r( 2 :; y) dz dy 


g( z ) ^ 2 k f( z ) ~ d{z) ~r- || k a 


1 


g 

1 - 


dz , 


where H-j^ represents the usual L 1 norm on I. The first condition of the theorem (Cl) guarantees that 


q(z) + -kf(z ) — p,(z) > 0 for all 2 £ / , 


so we can choose r in K r sufficiently small such that <&[/] > 0, i.e., $ : K r —>• K. On the other hand, 
using the assumptions (A1)-(A6), it is straightforward to show that 4>(A' r ) C C 1 (/)[^] This in turn, from 
Arzela-Ascoli theorem, implies that the operator $ has relatively compact image. 

Next we prove the second assumption of the theorem. For the sake of contradiction, suppose that there 
exist / £ K r with ||/|| M = r and A > 1 such that 


*[/] = A /■ 


2 Recall that continuous differentiability implies uniform continuity, and thus equicontinuity 
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Then it follows that 


A/p) < 


q(y) 

9 (y) 


f(y) dy - 


My)/ 2 + My) 


y(y) 


f{y ) dy 


f 0 J z 


r(g; y)kf(y) 

g(y ) 


f(y) dydz 


rx 1 

J X J z 


r(z; y)kf(y) 

g{y) 


t(y)dydz + \j o l M(p)dydz 


< 


J g{' : y ) J Q r ^’ ^ + dy 


g{y) 


g(y) + M/(y) - My) /(y) dy + ll/IL ■ IM 


< 


/ My) (y kf ^ L r ^ ; ^ + 9( ' 2/ ' ) ) dy 


My) 


g(y) + 2 fc /(y) - My)) dy + r ■ \\K 


1 

2 “ 

g 

1 


which yields that 
A < sup 


\\K 


My) 

1 

g 


My) J T (z-, y) dz + q{y)j dy + J 


g(y) 


My) + M/(y) - My) ) dy 


From the second condition of the theorem (C2) it follows that 

c Xi i 


sup 


y(y) 


M(y) J T (z]y)dz + q(y)^jdy + J (q{y) + \kf{y) - dy 


< 1. 


Then, we can choose r sufficiently small such that it contradicts the first assumption of the fixed point 
theorem, A > 1. 

Next we will derive conditions for the second condition of the fixed point theorem. For the sake of 
contradiction, let us choose k = 1 6 A'\{0} and assume that there exists / £ K r with ||/|| u = y < r and 
A > 0 such that 


/-$[/] = Xk = A. 


This equation in turn can be written as 


Afc = A = f(x) - $[/]p) 

f 1 My) 


< 


g(y) 


f(y) dy 


My )/ 2 + My) 


< rj — 


r P; y)kf(y) 

g(y) 

f(y) dy 4 


/(y) dydz + 


g(y) 

r m 


f(y )dy 
f Xl k a (z, y) 


M~) 


My) 

y(y) 


< 


My) ~ h f (y)/2 - q{y) 

My) 


1 My) - My )/ 2 

y(y) 

/(y) dy + q 2 ■ || k a 


g{y) 

f(y)dy + y 2 ■ II k, 

2 


/(y) dy dz 


V- 


(2.5) 


which should hold for all x £ I. Thus, provided that the condition (Cl) holds, we can choose rj £ (0, r) 
sufficiently small such that we get a contradiction to the second assumption of the fixed point theorem, 
A > 0. Hence, the fixed point theorem guarantees the existence of a positive fixed point of $ satisfying the 
bounds (|2.3|. 


Therefore, the function p* = /*/p is a stationary solution of the flocculation equations (1.11. Moreover, 
from the assumption (Al) and the continuity of the fixed point /* it follows that p* is non-zero, non-negative 
and continuous on /. □ 
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3. Principle of linearized stability and regularity properties of the linearized semigroup. 

In this section we summarize the principle of linearized stability as it applies to semigroups in general and 
our flocculation equation in particular. 

For a given autonomous ordinary differential equation, 


u = f(u ), 


the method for determining the local asymptotic behavior of a stationary solution it*, /(u*) = 0 , by the 
eigenvalues of the Jacobian Jy(n*) is quite well-known. In semigroup theory this method is known as the 
principle of linearized stability and was developed in the context of semilinear partial differential equations 
in 1331521155] ■ Later, Kato 137] extended this principle to a broader range of nonlinear evolution equations. 
Before presenting the principle of linearized stability we introduce some terminology, which can be found in 
many functional analysis books (see m for instance). 

The growth bound wo(A) of a strongly continuous semigroup (S(t)) t>0 with an infinitesimal generator 
A is defined as 


wo (A) := inf < w G 


3M U > 1 such that 
||5(f)|| < M u e ut for all t > 0 


The operator DA(f) denotes the Frechet derivative of an operator A evaluated at /, which is defined as 


DA(u)h, = A[u + h\ — A[u] + o(h), Vw G 2?(A), 

where o is little-o operator satisfying ||o(/i)|| < b(r) ||/i|| with increasing continuous function b : [ 0 , oo) —► 
[ 0 , oo), 6 ( 0 ) = 0 . 

The discrete spectrum ctd(A) of an arbitrary operator A on a Banach space A', is the subset of the point 
spectrum of A, 


f 7 p(A) = {A G C | 30 0 G X s.t. Acp = X(f>} 


such that A G <jd( A) is an isolated eigenvalue of finite multiplicity, i.e., the dimension of the set 

{ip G X : Aip = A ip} 

is finite and nonzero. Let (T (£)) t>0 be a Co semigroup on the Banach space X with its infintesimal gener¬ 
ator A. Then the limit wi(A) = lim^oo f -1 log (a[T(t)]) is well-defined and called the a-growth bound of 
(T(t)) f>0 . The function a[T(t)] is a measure of non-compactness of the semigroup T(t) as defined as in |38| . 
This measure associates non-negative numbers to operators (or sets), which tells how close an operator (or a 
set) is to a compact operator (or set). For example, for a bounded set M in a Banach space, a[M] = 0 implies 
that M (closure of M) is a compact set. Analogously, for a semigroup (T(f)) t>0 , a[T(t)\ = 0 indicates that 
the semigroup is eventually compact. 

With the above definitions, we are now ready to present the principle of linearized stability in the form 
of the following proposition (see [55] for the complete discussion of the proof of the following proposition). 

PROPOSITION 3.1. Define the nonlinear operator Af : D(F) C L 1 (I) —> L 1 (/) and let /* G D(M) be 
a stationary solution of (1-1), i.e., A /"[/*] =0. If Af is continuously Frechet differentiable on L X (J) and the 
linearized operator C = DAI (ft,) is the infinitesimal generator of a Cy-semigroup T(t), then the following 
statements hold: 

1. If wo ( C ) < 0, then /* is locally asymptotically stable in the following sense: There exists 77 , C > 1, 

and a > 0 such that if \\f — /*|| < 77 , then a unique mild solution T(t)f, satisfies\\T(t)f — || < 

Ce~ at ||/ — /*|| for all t > 0. 

2. If there exists Aq G cx(jC) such that Re A > 0 and 


max< wi(£), sup ReX><Re\o, (3-1) 

Ae<T£)(£)\{Ao} J 

then ff is an unstable equilibrium in the sense that there exists e > 0 and sequence {f n } in X such 
that f n /* and ||T(n)/„ - /*|| > £ for n = 1,2,... . 
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Having the explicit statement of the principle of linearized stability in hand, we now show that the 


nonlinear operator T defined in (1.4) satisfies all the conditions of Proposition 3.1 Towards this end, we 


first establish the elementary assumption of Propositio n |3.1| in the following lemma. 

Lemma 3.2. The nonlinear operator T defined in (1.4) is continuously Frechet differentiable on L 1 (/). 


Proof. The Frechet derivative of the nonlinear operator T is given explicitly as 

DF(<t>)[h{x)\ = -d x \gh](x) - (p(x) + ^k f (x)J h(x) + J T(x; y)k f (y)h{y) dy 
1 f x 

+ 2 J k *( x ~ y> v ) my)K x -y) + h(y)<t>i x - y )] dy 

nX\—X pXi—X 

~h(x) k a (x,y)<j>(y)dy-4>(x) k a (x, y)h{y)dy. 

J 0 J 0 

For the arbitrary functions ui, U 2 G L 1 (/) we have 

| DF( Ul )h(x) - DF(u 2 )h(x)\ < * Klloo |«i (y) - u 2 {y)\\h(x -y)\dy 


II K 


\h(x)\ || k a 


\h(y)\\ u i(x -y)- u 2 (x - y )| dy 


Wi{y) - u 2 (y)\dy 
[ \h(y)\dy 


+ KO) ~U 2 (X)\ | | kg, | | QQ 

Jo 

Consequently, taking the integral of both sides with respect to x and an application of Young’s inequality 
for convolutions (see |D Theorem 2.24]) to the first two integrals yields 

\\DT(ui)h(x) - DF(u 2 )h(x)\\ 1 < ||fc a |loo IK ~ Kll INIl + Halloo IK “ Kli IHli + IK “ Kll II^Hoo Mi 
for all h £ L 1 (/). Then it follows that 

II DF( Ul ) - DFMW, < 3 IlfealL K - Kli , 

which in turn implies that the nonlinear operator T is continuously Frechet differentiable on L 1 (/). □ 


In the previous section we have shown that the nonlinear operator T (1.4) has at least one non-trivial 


stationary solution, p* (in addition to trivial zero stationary solution). To derive stability results for this 
stationary solutions we first linearize the equation (1.1) around p*. A simple calculation yields that the 


Frechet derivative of the nonlinear operator T evaluated at a stationary solution p* (see Theorem 2.2) is 
given explicitly by 


C[h\(x) = DF(p*)[h\(x) = -d x (g{x)h{x)) - A(x)h(x) + f T(x; y)kf(y)h(y) dy 

J X 

—x r>x 

E(x,y)h(y)dy+ E(x - y, y)h{y) dy , (3.2) 

JO 

where 

E{x , y) = k a (x, y)p*{x) 

and 

^ rXx—X 

A{x) = -k f (x) + n(x) + J E(y , x) dy . 

We first prove that the linear operator C is an infinitesimal generator of a strongly continuous semigroup 
ST = (T(f)) t>0 . Consequently, we will prove two regularity results for the semigroup ST , which will prove 









useful in the spectral analysis of the operator £. Particularly, we will show that under some conditions on the 
model ingredients the semigroup ST is positive and eventually compact. The main implication of eventual 
compactness is that the Spectral Mapping Theorem holds (see |22| ) for the semigroup ST, 


<t (T(t)) \{0} = exp (ta(C )), t>0. 


Consequently, we will use the positivity of the semigroup ST in Section 5.2 where we employ the positive 
perturbation method introduced in [[25] , 

Lemma 3.3. If we define the domain of the linearized operator £ as 


V(C) = {</>€ L\I) | (#)' £ L\I), (#)(0) = m} , 

then the operator £ generates a Cq semigroup on 2?(£). 

Proof. The linear operator £ can be written as the sum of an unbounded operator 

Ci[h](x) = —d x (g(x)h(x)) — A(x)h(x) 


(3.3) 


(3.4) 


and bounded operators 


C 2 [h\(x)= T(x; y)k f (y)h(y)dy- 


E (x, y)h{y) dy, C 3 [h](x)= E(x - y, y)h{y) dy . (3.5) 


From the fact that g(x), A(x) £ C(I) and from the Lemma 2.4 of |5| it follows that C\ generates a Cq 
semigroup on 2?(£). Consequently, the bounded perturbation theorem of [46| §3, Theorem 1.1] yields that 
the operator C is also an infinitesimal generator of a Cq semigroup.□ 

Lemma 3.4. For a given stationary solution p * £ C(I) the operators C 2 : 2?(£) — ► L 1 (J) and £3 : 
2?(£) —> £*(/) defined in (3.5) are compact operators. 

Proof. We first prove that the operator £2 is compact. Then compactness of the operator £3 follows from 
analogous arguments. Let us denote a unit ball centered at zero in £ 1 (/) by B = {</> £ L 1 (I) | < l}. 

Recall that an operator is compact if it maps a unit ball into a relatively compact set. Consequently, observe 
that the assumptions (A2) and (A 6 ) together imply that the operator 


pXi —X 

d x C 2 {h\(x) = k a (x, xi - x)p„(x)h{x 1 - x) - d x (k a {x, y)p*{x))h(y) 

Jo 


dy 


+ 


d x T(x-, y)k f (y)h(y)dy-T(x; x)k f (x)h(x) 


is also bounded. Hence £ 2 [H] C H /1,1 (/) and from the Rellich-Kondrachov embedding theorem (see Q 
Theorem 6.3] for a statement of the theorem) it follows that the set £ 2 ]^] is relatively compact.□ 

Lemma 3.5. The operator £ defined in {3.tty generates an eventually compact Cq semigroup. And thus, 
the spectrum of the operator £ consists of isolated eigenvalues of a finite multiplicity only, i.e., cr(£) = <td(£). 


Proof. The operator £1 defined in (3.41 is well-known operator in size-structured dynamics literature. 
If g £ C 1 (/) and A £ C(I), then in Farkas and Hagen [26] Theorem 3.1] it has been shown that the Cq 
semigroup generated by the operator C\ is compact for t > 2 f Q 1 dy. The condition g £ C 1 (/) follows 
from our main assumption (Al), and continuity of the function 


1 !■Xl-X 

A (x) = ~jkf(x) + p(x) + J k a (x, y)p*{y ) dy 


follows from the assumptions (Al)- (A 6 ). Thus the semigroup generated by C\ is eventually compact. 
Conversely, in Lemma (3.41 we have shown that the operators £2 and £3 are compact. Hence, the Co 
semigroup generated by the operator £ = £1 + £2 + £3 is also compact for t > 2 f Q 1 dy. 

Therefore, the eventual compactness of the semigroup ST (generated by £) combined with Theorem 
3.3 of [ Til . §2.3] and Corollary 1.19 of [221 §4] together imply that the spectrum of £ consists of isolated 
eigenvalues of finite multiplicity.D 

Lemma 3.6. For a given steady state solution p* let us choose the functions k a , kf and T such that 
d x ( k a (x , y)p*(x)) < 0 for all x £ I and y £ ( 0 , x) 
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(3.6) 






and 


T(x; y)kf(y) > k a (x, y)p*{x ) for all x £ I and y £ [x, Xi) 


(3.7) 


Then the operator C generates a positive Cq semigroup. 

Proof. In Farkas and Hagen [26} Theorem 3.3] it has been shown that the operator C\ : £>(£) —► £ 1 (7) 
generates a positive Co semigroup under the main assumptions (A1)-(A6). On the other hand, from the 
conditions (3.6 1 and (|3.7|) it follows that 


C 2 [h]{x) + C 3 [h](x) > / k a (x-y,y)p*(x-y)h(y)dy- 
J o 

rx i 

- / k a (x, y)p*{x)h(y) dy 
Jo 


T(x; y)k f {y)h{y)dy 


> 


[k a {x - y, y)p*{x - y) - k a {x, y)p*{x)} h(y) dy 


+ / [r(z; y)kf(y) - k a (x, y)p*{x)] h{y) dy > 0, 

J X 

which in turn ensures that the operator C 2 + £3 is a positive operator. Since the positivity of a semigroup is 
invariant under a bounded and positive perturbation of its generator (see [22l §6, Corollary 1.11]), the result 
follows immediately. □ 


Remark 1. Lemma 3.6 has very important consequence. Specifically, if the positivity conditions (3.6 1 


and (3.7 ) hold and the spectral bound s(£) = sup {Re A | A £ cr(£)} is not equal to — 00 , then s(£) belongs to 
the spectrum <r(£) Engel and Nagel [ B21 §Theorem 1.10]. Moreover, the positivity and eventual compactness 
of the semigroup 3? together imply that the spectral bound s(£) is one of the eigenvalues of £ with finite 
multiplicity. 

4. Linearized stability and instability criteria for the zero stationary solution. In this section 
we will derive linearized stability results for the zero stationary solution of the flocculation equation. In 
contrast to non-trivial stationary solutions, zero stationary solution always exists (provided that the well- 
posedness assumptions (A1)-(A6) hold true). As we have discussed in Sectio n [ 3 ] t he stability of the steady 
states depends on the spectral properties of the linear operator £ defined in (3.21. We define the operator 
j\4 as the linear operator £ evaluated at the trivial stationary solution, p* = 0 


M[h]{x) = -d x [gh\(x) - 


MO) + 7jkf{x) 


h(x) + / T{x; y)k f (y)h(y)dy. 


(4.1) 


The assumptions (A1)-(A6) also ensure that the regularity conditions of Section [ 3 ] are all satisfied. 
Hence, the operator A4 generates a positive, eventually compact and strongly continuous semigroup. By 
Remark [l] we know that the spectral bound s(M) of the operator A4 is a dominant eigenvalue of Ai with 
finite multiplicity. Then, by the principle of linearized stability (Proposition ED* the stability of the zero 
stationary solution depends on the sign of this dominant eigenvalue. Thus, in the subsequent two subsections 
we derive conditions which guarantee positivity and negativity of the spectral bound s(M), respectively. For 
a more thorough discussion of the stability of zero stationary solution we refer readers to ]-H]. 

4.1. Instability of the trivial stationary solution. The operator A4 can be written as the sum of 
an unbounded operator 


Mi[h\(x) = - d x [gh\(x ) - (g(x) + -k f {x)j h(x) 


and a bounded operator 


rx 1 

M 2 [h\(x)= T{x;y)kf{y)h{y)dy. 

J X 


In J2Bj, authors have shown that the operator A4i generates a positive, eventually compact semigroup. 
Moreover, authors have shown that the spectral bound of A4i is positive if 


gp) 

ffP) 


exp 


/i(s) + \kf{s) 


g(s) 


ds ) dx > 1. 


(4.2) 
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On the other hand, we note that Ai 2 is a positive operator. Then, Corollary 1.11 of |221 §6] yields that 
the operator Ai = Aii + Ai 2 also generates a positive, eventually compact semigroup. Furthermore, the 
following inequality holds for spectral bound of A4i and Ai, 

s{A4 1 )<s(A4i+A4 2 ) = s{A4). (4.3) 


Consequently, this implies that the operator Ai also has a positive spectral bound provided that the condition 
(4.2) is satisfied. At this point, in Proposition 3.1 choosing Ao equal to the eigenvalue of Ai corresponding 
to s(Ai) and using Lemma 3.5 yields 


max < w i(Ai), 


sup 

ct £> (.M) \ { Ao } 


Re A > = 


sup 
AG<td(7W)\{A 0 } 


Re A < Re Aq . 


Then, the operator Ai satisfies all the conditions of Proposition 3.1 and thus results of this section can be 
summarized in the form of the following condition. 

Condition 1. Assume that the assumptions (A1)-(A6) hold true. Moreover, assume that 



g(z) 

d(x) 


exp 



fi(s) + ^kf(s) 
9(s ) 



dx > 1, 


then the zero stationary solution of the flocculation equation is unstable. 

4.2. Stability of the trivial stationary solution. In this section we will prove that under certain 
condition on model parameters we can ensure that the spectral bound of Ai is strictly negative. Since the 
positivity arguments that we used in the previous section cannot guarantee negativity of s(Al), we use a 
direct approach to prove that growth bound of Ai is strictly negative, oJo(Ai) < 0. To achieve our goal we 
use the following version of the well-known Lumer-Philips theorem (see for instance |22l §2, Corollary 3.6] 
and [T3] Theorem 2.22]). 

Theorem 4.1. (Lumer-Philips) Let a linear operator A on a Banach space (X, ||-||) the following are 
equivalent: 

1. A is closed, densely defined. Furthermore, A — XI is surjective for some A > 0 (and hence for all 
X > 0) and there exists a real number w such that A — u>I is dissipative, i.e., 

||/ - X(A - w/)/|| > ||/|| for all A > 0 and f £ V(A). 

2. Then, (A, T>(A)) generates a strongly continuous quasicontractive semigroup (T(t)) t>0 satisfying 


\\T(t)\\<e ut for i > 0. 


In the following lemma, we show the operator satisfies the first part of the Lumer-Philips theorem. 
Particularly, we establish that there exist a strictly negative real number w < 0 such that the operator 
M — ujI \s dissipative. 

Lemma 4.2. Assume that the assumptions (A1)-(A6) hold true. Then the linear operator Ai defined 
in (j.l) is closed, densely defined operator on the Banach space L 1 (I), and for sufficiently large A > 0 the 
operator Ai — XI : V(Ai) i—> L 1 (/) is surjective. Furthermore, if 


p.(x)-q(x) - i k f (x ) > 0 


(4.4) 


for all x £ I, then there exists a > 0 such that the semigroup (T(t)) t>0 generated by Ai satisfies the estimate 

||T(t)||i < e~ at for allt> 0. 


Proof. Since the operator Ai generates a strongly continuous semigroup (Lemma |3.3[ ), the first argument 
of the lemma is an immediate consequence of the Generation Theorem of |221 §2.3]. We now prove that there 
exist a > 0 such that Ai + ai is dissipative. For a given / £ V(Ai) = T)(C) and some h £ H and A > 0 we 
have 


f-X(Ai+aI)f = h. 


li 









Consequently, multiplying both sides by the sign function of / yields 


1/0)1 = /0) sgn (/(x)) 

= ~X [p(x)/(x)]' sgn (/(x)) dx- X 


-/c/(x) +/i(x) 


/ 0) sgn (/0)) 

+ A sgn (/0)) / r 0; y)kf{y)f{y)dydx + Xaf(x) + h(x) sgn(/(x)) , 


(4.5) 


where function sgn(/(x)) is defined as usual with sgn(O) = 0. For a given / £ V(A4) the set of points for 
which / does not vanish can be written as a finite union of disjoint open sets Ij = ( a,j , bj), i.e., /(x) yf 0 
for all x £ Uj =l Ij = (0, Xi). |^] On each interval Ij the function / can be either strictly positive or strictly 
negative. Moreover, on the boundaries we have f{a 3 ) = 0 and f(bj) = 0 unless a 3 = 0 or bj = X\. Then, 
integrating both sides of (|4.5[) on a given interval Ij = ( a,j , bj) we have 


f 


1/0)1 dx < -Xgibj) |/0j)| + Xg(aj)\f(aj)\ - X 


+ X 


"J 

r b o r x 

Ian Jx 


ILj L 

rbn 


-a+ -k f 0) +Mx) 


1/0)1 dx 


r(x; y)kf(y)\f(y)\ dy - 


Consequently, by summing (4.6) for j = 1,..., n we get 

r xi 


/0 


1/0)1 dx < -Ag(xi) |/0i)l + Aff(O) |/(0)| - A 

nX\ pXi 


\h(x)\ dx. 


~ a + 2 fc /0) + y{x) 


(4.6) 


1/0)1 dx 


+ X 
< -A 


'0 J x 

rx i 


r(x; y)kf{y) \f(y)\ dydx+ / \h(x)\ dx 


1. 


-a - q(x) + -fc/(x) +/i(x) 


1/0)1 dx 

+ A / 0(y)l/(y)l/ r 0 ',y)dxdy+l \h(x)\dx 

Jo Jo Jo 


-1 


= -A 

Hence, provided that we have 

for all x £ I, it follows that 


-a - q(x) - -fc/(x) + m0) 


rx i 

|/(x)| dx + / |/i(x)| dx . 

Jo 


-a - g(x) - -kf (x) + mO) > 0 


(4.7) 


li < 


= 11/ — A(A4 + al)f\\j 


In fact, if (4.41 holds true, then there exists a > 0 such that A4 + al is dissipative. Consequently, the result 


follows immediately from the Lumer-Philips theorem. □ As a direct consequence of Proposition |3.1| and the 
above lemma, we summarize the results of this section in form of the following condition. 

Condition 2. Assume that the assumptions (A1)-(A6) hold true. Moreover, assume that 

q(x) + \k f {x) - n(x) < 0 

for all x £ I, then the zero stationary solution of the flocculation equation is locally exponentially stable. 

5. Linearized instability and stability criteria for non-trivial steady states. In this section we 
present linearized stability results for the non-trivial stationary solution p* 0. We first derive conditions 
for instability (Section |5.1[) and then derive conditions for linear stability (Section |5.2[). 


3 See also M and m for similar partitioning in dissipativity proofs 
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5.1. Linearized instability. Recall that, from Proposition |3.1| instability of the non-trivial stationary 
solution depends on the spectral properties of the operator £. Specifically, the spectrum of £ contains at 
least one point Aq £ er(£) satisfying the instability condition (3.1). Towards this end (as we did in Section 


4.1), we first show that the operator £ has a positive spectral radius. 


Lemma 5.1. Assume that the positivity conditions (3.6)-(3.7) hold. Moreover, if the model parameters 
satisfy the following condition 


djx) 

9(x) 


exp - 


M s ) + \k f (s) + f 0 Xl s k a (s, y)p*(y) dy 


rx 1 — s 

Jo _ 

g(s) 


ds dx > 1. 


(5.1) 


then the operator £ has a positive spectral radius. 

Proof. Recall that the operator £ can be written as the sum of the operators £i, £2 and £3. Moreover, in 
Lemma [3.6|we have shown that the operator £1 generates a positive semigroup and the positivity assumptions 
(3.6) and (3.71 ensure the positivity of the operator £2 + £3. Therefore, from Corollary 1.11 of |221 §6] it 


follows that the spectral radius of £ is always greater than the spectral radius of the operator £3, i.e., 


s(£i) < s(£i + £2 + £3) — s(£). 

Conversely, provided that the condition holds, the arguments of j! 
that the spectral radius of the operator £1 is strictly positive. This result, combined with the inequality 


(5.2) 

Theorem 5.1] can be used to show 


(5.2) implies that the spectral radius of £ is strictly positive. □ We are now ready to present the main result 


of this section in the form of the following condition. 

Condition 3. Under the main assumptions (A1)-(A6) and the positivity conditions (3.6)-(3.7) the 


non-trivial steady state solution of the nonlinear evolution equation defined in (1.1) is unstable if 


gfo) 

g(x) 


exp - 


p(s) + \k f (s) + f Xl s k a (s, y)p*{y ) dy 


rxi — s 

Jo _ 

g(s) 


ds dx > 1. 


(5.3) 


Proof. Recall that from the proof of Lemma |5.1| it follows that 

s(£i) < s(£), 


where the operators £1 and £ are defined in (3.2). Note that Farkas and Hagen 
operator £1 has a positive spectral radius provided that the condition (5.3) holds true. 


have shown that the 
Consequently, if 

the condition (5.31 holds true, it follows that the operator £ has a positive spectral radius. Then from 
Proposition 3.5 and Remark [l] it follows that s(£) £ <td(£). Moreover, Proposition 3.5 together with [551 
Remark 4.8] imply that a-growth bound of £ is equal to negative infinity, oj\(C) = — 00 . Therefore, in 


Proposition 3.1 choosing Aq equal to the eigenvalue corresponding to s(£) yields 


max < wi(£), sup Re A > = sup Re A < Re Ao , 
I AGct£,(£)\{Ao} I AGcto(£)\{Ao} 


and implies that the non-trivial stationary solution p* of the nonlinear evolution equation defined in (1.1) is 
unstable. □ 

Remark 2. Let S C L X (I) be the set of non-trivial stationary solutions and S 1 (subset of S) denote 
the set of non-trivial stationary solutions existence of which guaranteed by Theorem \2.2\ For a stationary 
solution p * £ Si the modeling terms need to satisfy the conditions (Cl) and (C2). Consequently, plugging 
x = 0 into (C 2) yields the inequality 

f Xl q(x) 


Conversely, the instability condition 


J 0 g(x) 

implies that 


dx < 1. 


1 < 


g{x) 

g{x) 


exp 


y{s) + lk f (s) +J 0 X1 s k a (s, y)p*(y) dy 


rxi — s 

Jo _ 

g(s) 


ds dx < 


g(x) 

g{x) 


dx. 


which contradicts the existence condition (C 2). This in turn implies that stationary solutions in the set Si 
do not satisfy the instability condition. However, we note that Si is only subset of S, and thus the results of 
this subsection is only valid for non-trivial stationary solutions in the set S\Si- 
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5.2. Linearized stability. In Section |5TT| we have shown that the spectrum of the operator £ is not 
empty. This result, together with Proposition |3.6| and Remark |TJ imply that the spectral radius of £ is one 
of the eigenvalues of the operator £, so it is sufficient to show that all the eigenvalues of £ have a negative 
real part. However, to the best of our knowledge, the eigenvalue problem 

C[(j)\ = A <j) 


does not have an explicit solution. This forces us to utilize the positive perturbation method Farkas and 
Hinow |25j to locate the dominant eigenvalue of £. This method relies on the fact that compact perturbations 
do not change the essential spectrum of a semigroup. Towards this end we will perturb the operator £ = 
C\ + £2 +£3 (the operators £1, £2 and £3 are defined in (3.41-(3.5)) by a positive compact operator so that 


we can identify the point spectrum of the resulting operator. 

Lemma 5.2. Let us define the operator C as 

C[f](x) = d [ f(y) dy, 

Jo 

where ci = || k a • p*|L + ||T • fcylL- Then the operator C — £2 — £3 is positive and compact. 

Proof. It is easy to see that C — £2 — £3 is a positive operator, i.e., 

c[f}{X) - £2[/](*) - £3[/ID > [ [||*o. -Ph.il - - y, y)p*{x - y)\ f(y)dy 

Jo 

+ r Ol r • MIoc ^ F D y) fc /(2/)] M dy> 0 V/€ (l\I)) + . 

J X 

Conversely, C is a bounded linear operator of rank one, hence it is compact. Then the compactness of 


C — £2 — £3 follows from compactness of the operators £2 and £3 (see Lemma 3.41. □ 

Now define the perturbed operator V as V := £ + C — £2 — £3 = £1 + C. Then the eigenvalue problem 
for the operator V reads as 


Xf-V\f\ = Xf-C 1 [f\-C[f] = 0. 
This equation can be solved implicitly as 

t>cr. 

1 „ Ci 


fix) = U\ 


T{ A, x)g(x) + U2 g{x)T{\, x) J 0 


T{ A, s ) ds , 


(5.4) 

(5.5) 


where 


and 


££ = 


q{y)f(y) dy, 


U, = 


f(y) dy 


T(X, x) = exp Q y ). 


Integrating the equation (5.51 on / yields one equation for solving for Ui and U 2. Moreover, multiplying the 
equation (5.51 by q(x) and integrating over the interval / we obtain the second equation for solving for U\ 

(5.6) 


and U 2 - Consequently, these two equations can be summarized in the following linear system, 


U\Aii{\) + U 2 (.A12(A) — 1) — 0 
EM.<421 (A) — 1) + U 2 A 22 {\) = 0 


where 


Mi(A) — f 
Jo 

A 2 i (A) = f 
Jo 


1 


dx, A 12 (X) = [ 
Jo 

rrf, s , X dx, A 22 ( A) = / 
o T( A, x)g(x) J o 


Cl 


0 T i\ x)g(x) 
X1 q{x) 


0 g(x)T(\, x) J 0 

X1 c iq (x) 


T{ A, s) ds dx . 


0 g{x)T{\ 


-—r f T{\,s)dsdx. 
, X ) J 0 
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If the eigenvalue problem (5.4) has a non-zero solution, then there is non-zero vector (U\, tfy) satisfying the 


linear system (5.6). On the other hand, if there is non-zero vector (Ify, tfy) satisfying the linear system (5.6 1 , 


then the eigenvalue problem has a non-zero solution. Hence A € C is an eigenvalue value of the operator V 
if and only if 


K( A) = det 


Alii(A) 
^21 (A) — 1 


^4i2(A) — 1 
A 22 (A) 


— Alll(A)yl22(A) — (1 — Hi2(A)) (1 — AfyfyA)) — 0. 


(5.7) 


In structured population dynamics the function K is often referred as a characteristic function of eigenvalues 
of an operator, and similar characteristic functions have been derived in |491 ;22, 25, [26] [28]. The main 
advantage of having the characteristic function K is the task of locating the dominant eigenvalue value of 
the operator C reduces to locating the roots of the function K. Hence, in the following lemma we show that 
under certain conditions on the model parameters all the roots of the characteristic function K lie in the left 
half of the complex plane. 

Lemma 5.3. Under the conditions 


Aia(O) < 1, ^ 21 ( 0 ) < 1 (5.8) 

and 

K(0) < 0 (5.9) 

the function K does not have any roots with non-negative real 'part. Furthermore, the function K has at least 
one negative real root. 

Proof. It is straightforward to see that 

Ail (A) = AI 12 (A) = A 21 (A) = A 22 (A) = 0 as A —>■ 00 , 


so 


lim I\ (A) = —1. 

A—>-oo 

Moreover, observe that for i = 1,2 and j = 1, 2 the functions A i;J : ffi. —> R + are non-increasing, i.e., 


d\A ij ( A) < 0 . 


Consequently, for A > 0 from (5.8 1 we have 


1 — Hi2(A) > 1 — ^ 12 ( 0 ) > 0 


and 


1 — H2i(A) > 1 — AfyfyO) > 0. 

Conversely, differentiating K(X) for A > 0 yields 

K' = A 11 A 22 + AnA '22 ~h A' 12 (1 — A 21 ) ~hA' 21 (1 — ^ 12 ) < 0 . 

>0 >0 

Thus the function I\ restricted to real numbers is non-increasing. This in turn together with the condition 
(5.9) implies that the function K does not have any positive real root. 

Now for the sake of a contradiction, assume that there is Ai = a — b\ £ C with a > 0 and b fy 0 such 

that 


AfyAfy =0. 


(5.10) 


G(x) 



1 

g(y) 


dy, 


Let us define 








then for A = a — bi we have 


. , , f Xl cos[ 6 G(:r)] . f Xl sin [ 6 G(a;)] 

A n (Ai) = / L { dx + i I L , dx. 


I o T(a, x)g{x) 


/o 7\a, aOste) 


This in turn implies that 


—-i4ii(a) < Re An(Ai) < Au(a). 


(5.11) 


Analogous arguments yields similar inequalities for A i 2 , A 2 i and A 22 . On the other hand, if (5.101 holds 


true then using (5.111 it follows that 

An(a)Ai2(a) = |An(Ai)Ai 2 (Ai)| = |(1 — Ai 2 (Ai)) (1 — A 2 i(Ai))| 

> |1 — Re -A 12 (Ar)| |1 — Re A 21 (A 1 )| 

> (1 — A 12 (a)) (1 — A 2 i(a)) > 0 
Since K(X) is non-increasing for A > 0, from ( |5.9[ ) we have 

K(a) < A'(0) < 0 . 


(5.12) 

(5.13) 


The equation (5.131 is equivalent to 

Au(a)A 12 (a) < (1 — A 12 (a)) (1 — A 21 (a)) , 


which obviously contradicts the equation (5.121. Hence, the function cannot have a complex root with a 
non-negative real part. 

To establish the last statement of the lemma, observe that the function 1 — Ai 2 (A) is continuous and 
non-decreasing with 

lim 1 — A12 (A) = —00 . 

A —>•—00 

Conversely, from the condition ( |5.8| ) we have 1 — Ai 2 (A) > 0 for A > 0. Thus by Intermediate Value Theorem 
there is Ao < 0 such that 1 — A 12 (Ao) = 0. Consequently, evaluating K at A = Ao yields 

K (Ao) = Aii(Ao)A 2 2 (A 0 ) — (1 — Ai 2 (A 0 )) (1 — A 2 i(Ao)) = An(Ao)A 22 (Ao) > 0. 

Hence, the function K has at least one negative real root, which completes the proof of the lemma. □ With 
the above lemma in hand, we can now state the main result of this subsection in the form of the following 
condition. 

Condition 4. Suppose that the conditions 

K{ 0) = A n (0)A 22 (0) - (1 - A 12 (0)) (1 - A 21 (0)) < 0, 


a 12 (o) = (||fc 0 -P.IL + l|r 


kAJ L‘w)l exp G 


g(s) 


dsdx < 1 


and 


A 2i (0) = f 1 ^ exp f- f X ds) dx 
Jo g{x) \ Jo g(s) ) 


< 1 


hold true. Then, the non-trivial steady state solution p * is linearly exponentially stable. 
Proof. Lemma |5. 3 1 implies that the operator £\ +C has negative spectral radius, 

s(£ 1 T C) <c 0. 


Conversely, from Engel and Nagel [321 § 6 , Corollary 1.11], Proposition 3.6 and Lemma 5.2 it follows that 
s{£) = s(Gi + £-2 + £ 3 ) ^ s(Gi + £ 2 + £3 + C — £ 2 — £ 3 ) = s{£ 1 + C) < 0. 
Consequently, from [22] § 6 , Theorem 1.15 ] it follows that 

lo 0 (£) = s(£) < 0 . 


Hence, Proposition |3.1| yields that the non-trivial steady state solution p* is linearly asymptotically stable. 

□ 
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6 . Illustration of the results. In this section, we illustrate the theoretical development of the paper 
by giving explicit examples. 

Example 1. First, observe that p{t, x) = 0 is always a stationary solution of the flocculation model. 
Therefore, we only have to worry about the stability conditions derived in Section ^ For the maximal floe 
size, removal and renewal rates we choose terms similar to that of [21 

p(x) = 1, q(x)=b(x + 1) 

where b has yet to be chosen. Moreover, we assume that growth and fragmentation rates are proportional 
to the size of the floe, 


g{x) = x + 1 and kf(x) = 2x . 


Since the stability of the zero stationary solution does not depend on the aggregation rate and the post¬ 
fragmentation density function, we only assume that the functions k a and T satisfy the main assumptions 


(A2) and (A5), respectively. Consequently, plugging in this values into the instability condition (4.21 yields 


g{%) 


exp 


p(s) 


g(s) 


'i k ' is) -d,\dx = b 


I exp J 1 ds ) dx = 6(1 — e ) > 1. 


Thus, provided that b > 2, the zero stationary solution of the flocculation model is unstable. 


Conversely, plugging in the above parameters in the stability condition (4.41 one obtains 


q( x ) + - kf[x) — p{x) = bx + b — 1<0. 


Therefore, provided that we have 0 < b < we can guarantee the local exponential stability of the zero 
stationary solution. 

Example 2. In contrast to the zero stationary solution, positive stationary solutions do not always exist. 
Based on the analysis of Section [2j if the model parameters satisfy the existence conditions (Cl) and (C2) 
of Section [2] then the flocculation model possesses at least one nontrivial stationary solution. Towards this 
end, we choose an exponential growth rate and a uniform probability distribution for the post-fragmentation 
density function, 


g{x) = exp(— gce), T(ie; y) = X{ °' y] ^ ^ . 

V 

For other model parameters we choose linear rates, 


q(x) = b(x + 1), kf(x) = cx , p{x) 


c 


2 


x. 


A straightforward computation shows that the above model parameters satisfy the first existence conditions 
(Cl). Furthermore, plugging in the above rates into the second existence condition (C2) yields 


( 6 . 1 ) 


Hence, provided that the constants a, b and c satisfy the inequality S there exists a positive integrable 
stationary solution p*. Before we carry out the stability analysis for this stationary solution p*, the above 
functions should also satisfy the regularity conditions of Section [3j However, the regularity conditions of 
Section [3] depend heavily on the explicit form of the positive stationary solution in k a (x, y)p if {x). Since 
our analysis does not provide exact form of the positive stationary solution, we further assume that the 
stationary solution satisfies p* (x) > 0 for all x £ I. Consequently, we tailor the aggregation kernel such that 


it satisfies the positivity condition (3.71 and thus we proceed by choosing the aggregation kernel as 


d(l-x)(l~y) 

p*(x)p*(p) 
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k a (x, y) 











where the constant d has to be chosen sufficiently small such that 


cminp*(a;) > d. 

x£l 


This form of the aggregation kernel indeed solves the problem and all the regularity conditions of Section [3] 
are all satisfied. 

We are now in a position to derive conditions for the stability (and instability) of the positive stationary 
solution p*(x). Towards thi s en d, plugging the parameters into the instability conditions stated in Condition 
[3] and using the inequality (6.11 yields 




m(s) + |Ms) + fo s Ms, y)p*(y)dy 

g(s) 


dx < J 


<M dx = b ^l 

g(x) 


- 1 - 2 ae a 


< 1 . 


This in turn implies that these parameters do not satisfy the instability condition. However, since our 
instability conditions are only sufficient conditions, this does not imply that the positive stationary solution 
p* is always stable. 

For the stability of the positive stationary solution, we plug in the above parameters to the stability 
conditions stated in Condition [4] and obtain 


. .. e a (a — 1) + 1 , . . , a + e a — 1 — 2ae Q 

Ai 2 (0) < ci-^- < 1, H_2i(0) < b - y 

a z a- 


< 1 . 


( 6 . 2 ) 


where 


Cl = ||fc a -^lloo + ll r • k f lloo = —-TT + c < 2c - 


Moreover, the last stability condition, K{ 0) < 0, yields 


N , , , (e“ — 1) (a — 2 + (2a 2 — 3a + 2)e a ) 

v4n 0 • A 22 0 < ci--^- \ - f—L 

a 4 


< 1 - ci 


e a (a — 1) + 1 


1 - b 


a + e a — 1 — 2 ae° 


< (1-Hi 2 (0))(1-H 21 (0)). 


Conversely, since ci < 2c this yields another inequality for the parameters a, b and c, 


o (e“ - 1) (a - 2 + (2a 2 - 3a + 2)e°) ^ o _e°(a-l) + l 

Zc . ^ ( JL Zc 


l-b- 


- 1 - 2ae a 


(6.3) 


Thus, provided that we choose the parameters a, b and c such that the inequalities (6.1)-(6.3) hold, we can 
guarantee local exponential stability of this positive stationary solution. For the convenience of the readers, 
in Figure [6J] we have illustrated the stability region for the parameters a, b and c. Observe that for larger 
values of the parameter a the stability region for the parameters b and c significantly shrinks. As the value 
of the parameter a increases the growth rate decreases. Thus, from biological point of view, in order to 
keep the non-trivial stationary solution stable the values of the parameters b and c, associated to removal, 
fragmentation and renewal rates, should also decrease. 


7. Concluding remarks. Our primary motivation in this paper is to investigate the ultimate behavior 
of solutions of a generalized size-structured flocculation model. The model accounts for a broad range of 
biological phenomena (necessary for survival of a community of microorganism in a suspension) including 
growth, aggregation, fragmentation, removal due to predation, and gravitational sedimentation. Moreover, 
the number of cells that erode from a floe and enter the single cell population is modeled with McKendrick- 
von Foerster type renewal boundary equation. Although it has been shown that the model has a unique 
positive solution, to the best of our knowledge, the large time behavior of those solutions has not been 
studied. 

Using a fixed point theorem we showed that under relatively weak restrictions, which balance removal, 
growth, fragmentation and renewal rates, the flocculation model possesses a non-trivial stationary solution 


18 














1.0 



Figure 6.1. Stability region for the parameters a, b and c. Positive a,b and c values lying inside the solid assure local 
exponential stability of the positive stationary solution. 


(in addition to the trivial stationary solution). We used the principle of linearized stability for nonlinear 
evolution equations to linearize the problem around the stationary solution. This allowed us to infer stability 
of the stationary solutions by the spectral properties of the linearized problem. We then used the rich 
theory developed for semigroups, to derive the stability and instability conditions for the zero stationary 
solution. To derive instability conditions for the non-trivial stationary solution, we employed an eigenvalue 
localization method based on the well-known Krein-Rutman theorem. Lastly, we used compactness and 
positivity arguments to derive conditions for local stability of the non-trivial stationary solutions. 

Lastly, even though we showed that the flocculation model has at least one non-trivial stationary solution, 
our analysis does not state these stationary solutions explicitly. The nonlinearity introduced to the model 
by Smoluchowski coagulation equations, makes the task of finding explicit stationary solutions challenging 
even for constant model parameters. On the other hand, when only Smoluchowski coagulation equation 
is considered, it has been shown that the model has closed form self-similar solutions for constant and 
additive aggregation kernels pm 153) . Perhaps, under some conditions on the initial distribution and model 
parameters, solutions of the flocculation model also converge to self-similar profiles. Hence, as a future 
research, we plan to further investigate self-similar solutions of the flocculation model. 
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